Abrikosov flux-lines in two-band superconductors with mixed dimensionality 
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We study vortex structure in a two-band superconductor, in which one band is baUistic and 
quasi-two-dimensional (2D), and the other is diffusive and three-dimensional (3D). A circular cell 
approximation of the vortex lattice within the quasiclassical theory of superconductivity is applied 
to a recently developed model appropriate for such a two-band system [Tanaka et al 2006 Phys. Rev. 
B 73 220501(R); Tanaka et al 2007 Phys. Rev. B 75 214512]. We assume that superconductivity in 
the 3D diffusive band is "weak", i.e., mostly induced, as is the case in MgB2. Hybridization with the 
"weak" 3D diffusive band has significant and intriguing infiuence on the electronic structure of the 
"strong" 2D ballistic band. In particular, the Coulomb repulsion and the diffusivity in the "weak" 
band enhance suppression of the order parameter and enlargement of the vortex core by magnetic 
field in the "strong" band, resulting in reduced critical temperature and field. Moreover, increased 
diffusivity in the "weak" band can result in an upward curvature of the upper critical field near 
the transition temperature. A particularly interesting feature found in our model is the appearance 
of additional bound states at the gap edge in the "strong" ballistic band, which are absent in the 
single-band case. Furthermore, coupling with the "weak" diffusive band leads to reduced band gaps 
and van Hove singularities of energy bands of the vortex lattice in the "strong" ballistic band. We 
find these intriguing features for parameter values appropriate for MgB2 . 



I. INTRODUCTION 

Vortices are complex objects even in conventional type- 
II superconductors, and vortex structure can reveal un- 
derlying physics of thepairing interaction. The 40 K 
superconductor MgB2 Q is the best material discovered 
so far for studying multiple-band superconductivity. It 
can be well described by a two-band model that con- 
sists of the "strong" quasi-2D a band (energy gap « 7.2 
meV) and the "weak" 3D tt band (energy gap « 2.3 
meV), and there is evidence of induced superconductiv- 
ity 0, H, 0, 0, [1| and strong impurity scattering in the tt 
band. By tunneling along the c axis, Eskildsen et al have 
probed the vortex core structure in the tt band, and have 
found that the local density of states (LDOS) is com- 
pletely flat as a function of energy at the vortex centre, 
with no signature of localized states [![. Moreover, the 
core size as measured by a decay length of the zero-bias 
LDOS was found to be much larger than expected from 
Hc2- The existence of two length scales in the vortex 
lattice has also been suggested by the /zSR measurement 
0. The thermal conductivity measurement [To'| has de- 
tected a rapid increase of delocalized quasiparticles for 
field well below Hc2- Unusually large vortex core size 
and the extended nature of quasiparticle motion have 
also been found in the two-band superconductor NbSe2 
[O, [H, in, H, E, El ■ There is the possibility of multi- 
band superconductivity in a number of other materials 
[13, [H, IllH mi, m, m, and understanding the elec- 
tronic structure in the vortex state of multiband super- 
conductors is of great interest. 

Theoretically, the vortex state in a two-band super- 
conductor has been studied assuming both bands to be 
in the clean and dirty [2^ limit. With superconduc- 



tivity (mostly) induced in one band, it has been found 
that there are two different length scales associated with 
the two bands, and that the averaged zero-bias LDOS 
increases rapidly as a function of field in the "weak" 
band. Neither of these models, however, are suitable 
for describing many MgB2 samples, in which the a and 
TT bands are in the ballistic and diffusive limit, respec- 
tively. We have recently formulated a unique two-band 
model, in which one band ("cr band") is ballistic and the 
other ("tt band") is diffusive, and there is little interband 
impurity scattering [H, i]. This picture is appropriate 
for describing MgB2, in which two-band superconduc- 
tivity is retained even in "dirty" samples. (See exten- 
sive references in ref. [i|.) Having MgBz in mind, in our 
model it is assumed that the a and tt bands are quasi-2D 
and 3D, respectively, and superconductivity is mostly in- 
duced in the tt band through the pairing interaction with 
the cr band. We have examined the electronic structure 
around an isolated vortex and have found that the zero- 
bias LDOS in the tt band can have a decay length much 
larger than in the a band. A particularly intriguing fea- 
ture that emerges in our model is the possible existence 
of bound states at the gap edge in the ballistic band, in 
addition to the well-known Caroli-de Gennes Matricon 
(CdM) bound states [1^, in the vortex core. 

In this work, we extend our two-band model to de- 
scribe the vortex lattice and investigate the effects of in- 
duced superconductivity and impurities in the "weak" 
diffusive band on the electronic properties in the mixed 
state. We have found the presence of highly delocalized 
quasiparticles for relatively weak field, in agreement with 
the experiment on MgB2 p, [l^l and the earlier theoreti- 
cal works [2^l25j. Furthermore, band gaps and van Hove 
singularities of energy bands of the vortex lattice in the 
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"strong" ballistic band can be reduced by coupling with 
the "weak" diffusive band. There exist gap-edge bound 
states in the vortex core in the ballistic band - that solely 
arise from coupling to the diffusive band - also in the vor- 
tex lattice. 



II. FORMULATION 

We utilize a recently developed model appropriate for 
a system with coupled ballistic and diffusive bands [H, [3] ■ 
Both the ballistic and diffusive limits can be described 
within the quasiclassical theory of superconductivity, in 
which all the physical information is contained in the 
quasiclassical Green function, or propagator. We write 
the propagator in band a as ga{c,PFa, R), where e is 
the quasiparticle energy measured from the chemical po- 
tential, pFa the quasiparticle momentum on the Fermi 
surface corresponding to band a, and R the spatial coor- 
dinate. The hat refers to the 2x2 matrix structure in the 
Nambu-Gor'kov particle-hole space. In the clean a band, 
(Jai^iPFa: R) Satisfies the Eilenberger equation p?!. [28j. 
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(1) 



where vp„ is the Fermi velocity and A^- the (spatially 
varying) order parameter. The three Pauli matrices are 
denoted by f^, i = 1,2,3, and [...,...] denotes the com- 
mutator. Motivated by the Fermi surface of MgB2, we 
assume a cylindrical Fermi surface and treat the a band 
as quasi-two-dimensional. The cr-band coherence length 
is defined as fo- = hvp^/2T:Tc, where Tc is the transition 
temperature, and used as length unit. 

We assume that the n band is in the diffusive limit. In 
the presence of strong impurity scattering, the momen- 
tum dependence of the quasiclassical Green function is 
averaged out, and the equation of motion for the result- 
ing propagator gT^{e,R) reduces to the Usadel equation 



era - A^r, (/tt 
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(2) 



Throughout this 
ij = D6ij and de- 



with the diffusion constant tensor B 
work, we assume an isotropic tensor 
fine the 7r-band coherence length — \/hD /2ttTc. Both 
ballistic and diffusive propagators are normalized accord- 
ing to gl=gl = -nH ^2% 

The quasiparticles in different bands are assumed to 
be coupled only through the pairing interaction. Selfcon- 
sistency for the spatially varying order parameters in the 
two bands is achieved through the coupled gap equations, 



A„(i?) = ^14/3^F/3^/3(i?), 
/3 



(3) 



where a,(3 £ {a, tt}, and A^ = fi Re Aq,— f2 Im A„. The 
coupling matrix Vaf3 determines the pairing interaction. 



Npf} is the Fermi-surface density of states on band /3, and 



^a{R) = / ^ (/<T(e,PFa,i?))p-^„tanh( ^ j, 
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fjr (e, R) tanh ( 



(4) 



Here fa is the upper off-diagonal (1,2) element of the 
matrix propagator ga, and Cc is a cutoff energy (that will 
be ultimately eliminated as discussed below). The Fermi 
surface average over the cr-band is denoted by (• • • )pp^. 
In this work, we consider isotropic s-wave coupling as is 
the case for MgB2. 

To study the mixed state, we introduce a circular cell 
approximation, meaning that we simulate the vortex unit 
cell by a circle. Other suggestions for circular cell approx- 
imations have been introduced in refs. [s^] and [3l| for 
the Usadel and Eilenberger equation, respectively. We 
assume strong type-II superconductivity, in which case 
the spatial variation of magnetic field within the unit 
cell can be neglected for fields not too close to Hd ■ This 
is the case for MgB2. We assume that magnetic field is 
in the — z direction (taking into account that the electron 
charge e < 0, this gives a positive phase winding). The 
vector potential is then given by 
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(5) 



with Hq = -HqCz = W X a and r = ^/x^^^^ . The 
radius of the vortex unit cell is determined by the fact 
that one flux quantum penetrates the unit cell. 



with $0 = /ic/2|e| 



7r?ic/|e 
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Thus, 
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(6) 



(7) 



This can be used directly in both the Eilenberger and the 
Usadel equation. In both cases, the vector potential can 
be incorporated by the replacement 



he 



T3A^,X 



(8) 



where the o symbol involves a time convolution if X is 
time-dependent; otherwise it simply implies matrix mul- 
tiplication. Note that in the diffusive case this also affects 
the expression for the current density. This can be ob- 
tained by the requirement of local gauge invariance that 
the local gauge transformation 
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dtX 



(9) 
(10) 

(11) 
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with any xi^^, t) should leave the basic equations of mo- 
tion (transport equations) invariant (in our case there 
is no time dependence, hence no scalar potential $). 
The vector and scalar potentials can be gauged locally 
away, by writing down the equations at any point R with 
Vx = "■^^ 9,nd dtx = f Thus, the equations with 
potentials can be obtained by observing that 



(12) 



Now, we use the Riccati parametrization of the Green 
functions for both the Eilenberger [s^, [s^, [sE] and 
Usadel j36|] equations: 



9a = 



1 - 7q7q 

1 + 7a 7a V 7„7„ - 1 y ' 



(13) 



where %{€,pFa,R) = j*{~€*,-pFcr,R) and 7^(e,i?) = 
7^(-e*,^)- 

Because the transport equations for the Riccati am- 
plitudes must also be invariant under any local gauge 
transformation, we can write them down immediately by 
noting that the above gauge transformation means 



7 — > e^^ 070 e*"^ 
7 —I- e~*-*^ o 7 o e"*-^. 



or if we have no time dependence. 
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(14) 
(15) 



(16) 
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Similarly, A ^ e^'^A and A* -> e^'^'^A*. We thus 
obtain for the clean a band, 

A, + 2e-f, + Aljl+ihvf,,,daa^0, (18) 
with di = (Vi + 2iVx) = (Vi - ijTcAt), and 

Al -2e% + A^ja+ifi^faA'f'y = (19) 



with di = (Vi - 2iVx) = (V, + ij^A,). For the dirty tt 
band, we have [391 



A^ + 2e-f^ + AIjI - ihD d,da^ 



and 



a; - 26 7^ -f a^jI + ihD dAi^ 



27,r(i9»7^)^ 

1 + IttItt 



2j^{di%)'^ 

1 + 7x7;r 



= 0, 
(20) 



(21) 

We consider a vortex extending in z direction, with its 
centre situated at a; = y = for each z. We write {x,y) = 
r(cos sin(/)) and (vFx,VFy) = fF(cos0, sin0), and A = 
|A(r)|e^'^, 7 = 7(6';0,r)e"^, and 7 = 7(6*; 0, r)e-*'^. 
Then we can use the above equations with the quan- 
tities |A(r)|, 'j{6; (f),r), and ^{9;(j),r), and with the re- 
placements 9 (V — ij^A^{r)e^ + ie^/r) and d 
[y + ij^A^{r)e^ — iSff,/ r) . Here we have used \/cf) — e^/r. 

To obtain the Riccati amplitudes in the 2D ballistic 
band, we must first solve for the boundary values along 
the unit cell circle, from which to integrate the Riccati 
equations along a given trajectory. To simulate the phase 
change across a unit cell boundary, we impose the bound- 
ary condition (here we omit the subscript cr). 



7i„(6'; (/),rc) = --7o«t(6'; </> + tt, r^) , 



(22) 



where in and out refer to incoming and outgoing trajecto- 
ries, respectively, and similarly for 7. For each trajectory 
we solve for the boundary value 7m selfconsistently. 
For the 3D diffusive band, we use the symmetry 



7^(0, r)e'''^ 
7^(0, r)e-'''^. 



(23) 
(24) 



Then equations (PO)) and ipT|) reduce to the following 
dimensionless equations for 7,^(0, r) — j^rir) = Jtt and 
7^(0, r) = %{r) = %: 



|A,|(l-f 7^)4- 26 7, -i (^^^ (d^^j^ + \j^ 



r 1 \ 1 - Jt,% 

— - - It:— ^ 

r 1-1-7^7^ 



277r(^r77r)^ 
1 + ItxItx 



0, 



(25) 



A^Kl-f 7^) -26 7^-1-1 



~ l-77r77r 27^(5r77r)^ 



1 + 7x7;r 



1 



= 0, 



(26) 



The boundary condition in the dirty case is that the first the (cutoff dependent) interactions and thus transform- 
derivative of 77r (r) and 7^ (r) is zero at the cell boundary. 



We solve for the order parameters by diagonalizing 
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ing the gap equations ^ into the form 

A(") \ f 
Ad) i I 



Po = 0.3 n^ = n„ T = 0.1T,: r^^„ = 5 



A(i) 



(27) 



where ^("^ and J-"^) are the anomalous amphtudes in 
the diagonal basis. The larger eigenvalue, say, A^^^, de- 
termines Tc and can be eliminated together with Cc- The 
smaller eigenvalue A^) can be parametrized by the cutofF- 
independent quantity 



A 



A(°)A(i) 
A(") - A(i) 



(28) 



We solve the set of equations (|18|19|25|26p along with the 
gap equations for the order parameters selfconsistently. 

After selfconsistency has been achieved, the LDOS in 
each band can be calculated by 

N„{e,R)/NFa = - Im {g„{e,pFa,R))pp„/T^, 
N^{e,R)/NF. = - Im5,(e,i?)/7r, (29) 

where ga is the upper diagonal (1,1) element of Qa- 

The current density around the vortex has contribu- 
tions from both the a and the tt band. The corresponding 



expressions are 




t{R) 


/•°° de 


2eNFa 


J-oo 27r 




D r°° 


2eNF^ 


I" i-oo 



{vfAt^ ga)pp„ tanh 



de 



\2TJ ' 



2^Im[/:9Mtanh(^— j. (30) 



In our model, the bulk behaviour of the system is com- 
pletely specified by four material parameters, po, n-^/ricr, 
Tc, and A. Here po is the zero-temperature bulk gap ra- 
tio, Po = l^bulk I ^bulk ^ ^ NFcKNFrr + NfjA. 

For MgBa, po ^ 0.3 and n^/ua- « 1 - 1.2 [1, 0, Eg- 
The zero-temperature gap equations for a homogeneous 
system relate A with po, n-^/na, and the bulk gap ratio 
near Tc pj. From this relation we find that for MgB2, 
A can be positive or negative [H, [s^ HO, |4l| and can be 
close to zero [s^. Negative A implies that the effective 
Coulomb interaction dominates over the effective pairing 
interaction in the subdominant A^^ channel, and in this 
case superconductivity is purely induced in the tt band 
(with > 0, while A^^) < 0). 

In the presence of inhomogeneity, there is another ma- 
terial parameter, namely i^tt/Cct, where ^o- — 6.8 nm 
(Tc — 39K) for MgB2 As discussed in the next sec- 
tion, to reproduce the experimental data of ref. Q we 
find ^7r/C(T ~ 2, for po = 0-3 and riT^jn^ = 1. The condi- 
tion for the TT band to be in the dirty limit so that the 
Usadel equation is applicable is ^jr/Co- < 5 for MgB2 0. 
We present results for po = 0.1, 0.3, 0.5; ^tt/^ct = 1, 3, 
5; nT^/n„ — 1, 1.2; and A = —0.1, 0, 0.1; and discuss 
the effects of induced superconductivity and hybridiza- 
tion of the diffusive and ballistic bands. For a given cell 
radius Tc, the corresponding field strength is given by 
equation Q in units of fi,c/|e|^^, which is about 13 T for 
MgBa. 




FIG. 1: Order parameters |A^(a;)| and |Act(s)| and their ratio 
as a function of coordinate x along a path through the vortex 
centre for ^,r/C<T = 1,3,5, po = 0.3, n-^ = n^, T = O.lTc, 
rc/ia = 5, for A = 0.1 (a,c) and -0.1 (b,d). 



III. RESULTS 
A. Order parameter 

In figure [T] we show the order parameter magnitudes 
in the two bands (a,b) and their ratio (c,d) as a function 
of coordinate x along a path through the vortex centre 
for Po = 0.3, tItt = Ua, and T = O.lTc, for various values 
of the coherence length ratio Ctt/Co-; for A = 0.1 (a,c) 
and -0.1 (b,d). The unit cell radius rc/£,a — 5 corre- 
sponds roughly to 0.5 T for MgB2. The vortex structure 
is affected significantly by the Coulomb interaction and 
the impurity scattering rate in the tt band. For A > 0, 
I A^ (a;) /Act (x) I is smaller than the bulk gap ratio po in 
the vortex core. While for ^,r/Co- = 1 the ratio recovers 
more or less to the bulk value at the cell boundary, as 
Ctt/^ct increases, |A7r(a;)| is suppressed more strongly and 
|A7r(a;)/ACT(a;)| is smaller than po in the entire unit cell 
(for ~ fc). Furthermore, through coupling with the 
dirty tt band, suppression of the (T-band order parameter 
by magnetic field is also enhanced for larger ^,r/Ccr- This 
suppression of the order parameter and enlargement of 
the core area in the a band with increasing S,Ti/^a are 
more drastic when the Coulomb repulsion is dominant in 
the TT band (A < 0). Interestingly, in this case, the deple- 
tion of |A7r(a;)| is relatively small and |A^(a;)/Acr(a;)| is 
substantially larger than po, especially for larger ^,r/Co-- 

Changes in the order parameter magnitudes in the two 
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^^/^a = 5 Po = 0.3 n^ = n„ T = 0.1T, 
2.0 



Po = 0.5 n^ = n„: r^5„=10 




(c) r /4„=6 




(d) r /5„=6 



1 2 3 4 



12 3 4 



FIG. 2: Order parameters jA.^(a;)| and |Act(x-)| and their ratio 
as a function of coordinate x along a path through the vortex 
centre for C,r/C<T = 5, po = 0.3, = n^, T = O.lTc, rc/C<T = 
4, 6, 10, for A = 0.1 (a,c) and -0.1 (b,d). 
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FIG. 3: Order parameter magnitude at the vortex unit cell 
boundary in the two bands as a function of temperature T for 
po = 0.3, Utt = n^, K = ±0.1, rc/^a = 10, for (a) = 3 

and (b) 5. The result for a single clean band is also shown. 
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FIG. 4: Same as figure |3] except for po = 0.5. 



bands and their ratio as the field strength changes are il- 
lustrated in figure [21 for ^tt/^o- = 5, po = 0.3, n-,^ — n^, 
and T = O.lTc, for various values of r^, for A = 0.1 
(a,c) and -0.1 (b,d). The result for rc/£,a = 10 is similar 
to that for an isolated vortex. It can be seen clearly 
that for A > 0, the 7r-band order parameter is sup- 
pressed strongly by magnetic field and the core area is 
enlarged. In contrast, for A < 0, the depairing effect is 
more manifest in the a band, with |A^(a;)/Acr(a;)| reach- 
ing about two times po in the vortex centre. It is clear 
that, while superconductivity is (mostly) induced in the 
TT band. Coulomb interactions can renormalise substan- 
tially the length scales and the core sizes in the two bands 
in different ways. 

We illustrate the influence of the diffusivity in the tt 
band on the critical temperature in figures [3] {po — 0.3) 
and [3] (po = 0.5) for = and rc/S,a — 10, for (a) 
^v/^tT = 3 and (b) 5. In these figures the order parame- 
ter magnitudes at the cell boundary, |A(r'c)|, in the two 
bands are plotted as a function of T for A = 0.1 and 
-0.1, along with that for a single clean band. Points are 
results obtained by selfconsistent calculation and curves 
are guides to the eye. In a single ballistic band, the criti- 
cal temperature where |A(rc)| vanishes is slightly smaller 
than the zero- field Tc- When the ballistic a band is cou- 
pled with the diffusive tt band, the critical temperature 
is reduced further, and superconductivity is more sup- 
pressed by magnetic field for A = —0.1 than for A = 0.1. 
This difference is enhanced for larger ^tt/Co-, as a result 
of stronger suppression of the cr-band order parameter 
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FIG. 5: Vortex core size as defined by equation (|3ip in the two 
bands as a function of temperature T for ^w/£,a = 3, po — 0.3, 
Un = n„ , A = 0.1, for (a) rc/^a = 10 and (b) rc/^o- ~ 5. The 
results for a single clean and dirty band are also shown. 



as discussed above. These effects of the diffusivity and 
the Coulomb interaction in the tt band can be drastic 
when the couphng between the two bands is strong, as 
demonstrated for po = 0.5 in figure [H 

A quantity to characterise the vortex core structure is 
the vortex core size defined by [42, 



.-1 ^ dAir 



0) 



1 



dr A{r 



(31) 



where r is the radial coordinate measured from the vor- 
tex centre, and A(r — oo) = A(rc) is the 'bulk' order 
parameter in the vortex lattice. Around an isolated vor- 
tex in a single clean band, the order parameter exhibits 
the KP effect 42], i.e., shrinkage of the vortex core size as 
T is lowered, approaching zero in the zero-temperature 
limit [il, 0, EBI . This is due to depopulation of higher- 
energy bound states in the vortex core. In an s-wave su- 
perconductor with nonmagnetic impurities, however, the 
core size as defined above saturates as T approaches zero 
[3, SB] • This stems from broadening of the bound core 
states that removes the singular behaviour in the spatial 
variation of the order parameter in the vortex core. The 
vortex core shrinking ceases when ksT becomes smaller 
than the energy width of the zero-energy bound states in 
the core. 

When a ballistic band and a diffusive band are coupled, 
the KP effect is induced in the diffusive band, as found in 
refs. [5, for an isolated vortex. In the vortex lattice, 
when vortices are well separated, one finds the KP effect 
in a single clean band as well as in coupled clean and 
dirty bands. This is demonstrated in figures [5] (A = 0.1) 
and [6] (A = —0.1), in which in the two bands is plotted 
as a function of T for ^-k/^o- = 3, po = 0.3, and = n^, 
for (a) Tc/^a — 10 and (b) rc/^a = 5- The core size in the 
case of a single clean and dirty band is also shown. As 
can be seen in figurc[5lja), for A > 0, the core size in the tt 
band is larger than that in the a band for all temperature. 



FIG. 6: Same as figure [5] except for A — —0.1. 
£/£ =3 p„ = 0.3 n, = n„ T/T =0.1 
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FIG. 7: Vortex core size as defined by equation in the 
two bands as a function of field for £,w/^,(y = 3, po = 0.3, 
= Ua, T/Tc = 0.1, for (a) A = 0.1 and A = -0.1. 



The TT-band order parameter exhibits the KP effect also 
when the Coulomb repulsion dominates (figure[Sl[a)), and 
in this case, as defined above is always smaller in the 
TT band than in the a band (except for T very close to the 
critical temperature). Also note that, with dominating 
Coulomb interactions, the T-linear behaviour of the KP 
effect is better developed in the tt band than in the a 
band. 

As magnetic field increases and vortices come closer 
together, the bound state wavefunctions of neighbouring 
vortices begin to overlap and form energy bands, and 



quasiparticles can travel through the periodic array of 
vortices [ill, IH) Ell (see also references in ref. [HI). 
We have performed our calculation for the cell radius 
Tc > 4. For this range of Tc, we find that in a single diffu- 
sive band, the core size becomes smaller as Tc decreases 
for any given temperature. Such shrinkage of the vortex 
core as field increases can be understood as due to inter- 
vortex transfer of quasiparticles jllj , with higher-energy 
core bound states turning into extended states. This is 
also a trend for a single clean band as well as coupled 
ballistic and diffusive bands for relatively high T and rel- 
atively large Tc- However, as can be seen in figures [SJb) 
and [6{b), for relatively strong field, is finite in the 
zero-temperature limit; since the derivative of the order 
parameter at the vortex centre remains finite as temper- 
ature approaches zero. Thus for a given (low) tempera- 
ture, can increase as a function of field strength above 
a certain critical value. This is illustrated in figure [3 
where is plotted as a function of field strength iJp for 
T/Tc = 0.1, ^7r/C<T = 3, po = 0.3, and = no-, for 
(a) A = 0.1 and (b) A = -0.1. In both bands, as field 
increases, the core size is reduced initially, but at some 
critical strength it starts increasing. In a single ballistic 
band, behaves similarly as a function of H^: this is 
consistent with the work by Miranovic et al [4§| , who ex- 
amined the effects of impurities on the vortex core size by 
including impurity scattering in the Eilenberger equation 
(single band). They have also found that for very small 
mean free path, the core size decreases monotonically as 
a function of field strength, as we find for a single dif- 
fusive band. In the case of two bands, in the tt band 
has nonmonotonic behaviour as in the a band, and the 
low-temperature value of is larger (smaller) in the it 
band than in the cr band for A > (< 0). 



B. Spectral properties of diffusive band 

In the diffusive tt band, the LDOS is completely fiat as 
a function of energy at the vortex centre, as found for an 
isolated vortex P, T^l and for coupled diffusive bands ^ | . 
This is in agreement with the observation with STM [5 1 . 
The vortex core size as measured by a decay length of the 
zero-bias LDOS has also been probed by STM When 
superconductivity is (mostly) induced in a "weak" band 
in a coupled two-band system, such a decay length can 
be much larger in the "weak" band than in the "strong" 
band [ij, H, HI, HH]. We find that in the vortex state, 
there is little change in the zero-bias LDOS as a function 
of r in the a band from that for an isolated vortex, ex- 
cept for rather small Vc- In contrast, as Tc is reduced, the 
TT band starts exhibiting substantial overlap of quasipar- 
ticle states across the cell boundary for relatively large 
Tc. Furthermore, the LDOS decay length in the tt band 
strongly depends on pQ and £,tt/£,<t- In figure [8] we present 
the zero-bias LDOS in the tt band as a function of dis- 
tance r from the vortex centre for tt-tt = jt-o-, T/Tc = 0.5, 
A = 0.1,0,-0.1, and rj^^ = 20; for (a) po = 0.3 for 
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FIG. 8: Zero-bias LDOS in the tt band as a function of dis- 
tance r from the vortex centre for n^r = ria, T/Tc = 0.5, 
A = 0.1, 0, -0.1, rc/e<T = 20, for (a) po = 0.3, C^/C<. = 1, 3, 5 
and (b) C,r/C<T = 3, po = 0.1, 0.3, 0.5. 
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FIG. 9: Half-decay length ^i, of the zero-bias LDOS in the tt 
band as a function of ^,r/^CT for po = 0.3, = no-, T/Tc = 0.1, 
Vc/^a = 20, A = 0.1, 0, -0.1. For MgB2, Cc2 = 10 nm and the 
observed ^„ ft; 30 nm corresponds to ^tt/C^ « 2 in our model. 



various values of ^tt/^ct and (b) ^tt/Co- = 3 for several val- 
ues of pq. For parameter values appropriate for MgB2 
(Ctt/^ct ~ 2, Po = 0.3; see below), the overlap of core 
states is not negligible already for rc/£,a = 20. For a 
fixed Pq, the diffusivity in the tt band (larger ^tt/Cct) en- 
hances intervortex transfer of 7r-band quasiparticles. For 



8 



5./^a = 3 n,/n„=1.2 TAT, = 0.5: r,/§„ = 20 




p = 0.3 
p = 0.5 
1-band cl?an 




(b)A = -0.1 






p = 
p^ 

1-band cj 


0.3 

0.5 - - 
3an 























E/T, 



FIG. 10: Vortex centre spectra in the a band as a function of 
energy e for po = 0.3, 0.5, ^n/£,cT — 3, rin/ncr = 1.2, T/Tc — 
0.5, rc/C<T = 20, for (a) A = 0.1 and (b) -0.1. The result for a 
single ballistic band is also shown. 




FIG. 11: Vortex centre spectra in the a band as a function 
of energy e for ^tt/^ct = 3, 5, = tIo-, A = —0.1, T/Tc = 0.5, 
rc/^iT = 10, for (a) po ~ 0.3 and (b) 0.5. The result for a 
single ballistic band is also shown. 



a given ^-w/^a-, the weaker the induced superconductivity, 
the stronger the extension and overlap of quasiparticle 
states. The Coulomb repulsion in the tt band (negative 
A) can enhance these effects further: e.g., for the param- 
eter set for figure [S^b), the dependence of the LDOS on 
A is significant for po = 0.1. 

To characterise the core size, we plot in figure IH] the 
half-decay length of the zero-bias LDOS in the tt band, 
^v, in units of ^c2 = 10 nm (the Ginzburg-Landau co- 
herence length in MgB2 Q) as a function of ^Tr/^a for 
Po = 0.3, — n„, T/Tc = 0.1, and r-c/^o- = 20. The 
^1, is a linear function of ^tt/£,(j for ^tt/Cct ^ 3 in this ex- 
ample, and increases rapidly as ^tt/Cct increases further. 
Also for larger ^jr/Co-, its dependence on A becomes no- 
ticeable. For MgB2, the experimental value of ^t, sa 30 
nm corresponds to ^tt/Co- ~ 2 in our model. In ref. [1^ 
such a plot of was made for coupled dirty bands for 
£,Tr/^a < 2.5: overall f„ in this case is slightly smaller 
than our ^y. 



C. Spectral properties of ballistic band 

A particularly interesting feature found in our model 
of coupled ballistic and diffusive bands is that, through 
coupling with the "weak" diffusive band, there can be 



additional bound states at the gap edge in the "strong" 
ballistic band. We find such bound states for an isolated 
vortex [U, as well as for the vortex lattice. In figure [TU] 
the LDOS at the vortex centre in the a band is plotted 
as a function of energy e for po = 0-3, 0.5, ■Ctt/'Cct = 3, 
n^/n„ = 1.2, T/Tc = 0.5, and r^/Ca = 20; for (a) A = 
0.1 and (b) -0.1. We present in figure [U] the LDOS for 
S,7r/Co- = 3, 5, n^r = n^, A = -0.1, T/T^ = 0.5, and 
rc/£.a = 10; for (a) po = 0.3 and (b) 0.5. The result 
for a single ballistic band is also shown in these figures. 
The peak at e = reflects the CdM bound states that 
arise from repeated Andreev scattering from the order 
parameter in the vortex core. In the core of an isolated 
vortex in a single clean band, the spectrum shows neither 
coherence peak nor additional bound state at the gap 
edge. In the vortex lattice, as quasiparticles travel in the 
periodic pairing potential, the continuum energy levels 
above the energy gap are grouped into bands and van 
Hove singularities appear at the band edges [48]. This 
can be seen clearly for a single ballistic band in figures fTOl 
andfm in the latter, the band widths are larger due to 
reduced lattice spacing. 

When the ballistic a band is coupled with the diffusive 
TT band, the vortex core area is enlarged and the peri- 
odic pairing potential becomes effectively weaker. As a 
result, band gaps and van Hove singularities are reduced 
as can be seen in figures [TU] and [TT] At the same time, 
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FIG. 12: Current density in the two bands as a function of 
distance r from the vortex centre for ^ir/Co- = 3, po = 0.5, 
n-n = ria, A = -0.1, rc/Ca = 10, for T/Tc = 0.1, 0.3, 0.5. 



extra bound states tend to appear at the gap edge. These 
features become prominent for larger po (0.3 vs. 0.5 in 
figure [TU)) and £,ttI£,<j (3 vs. 5 in figure [TT|l. and are more 
enhanced for negative A than for positive A (figure [TU)) . 
In most of the examples shown in these figures, the spec- 
trum above the energy gap is more or less continuous. 
Furthermore, in the two-band results presented in fig- 
ure [TT](b), the energy gap is reduced substantially and a 
band of gap edge bound states is formed. 
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FIG. 13: Magnetic field distribution 5H — H — Ho as a 
function of distance r from the vortex centre for ^-n/^a = 3, 
fin = A = -0.1, T/Tc = 0.1, rc/^a = 10, for po = 0.3, 0.5; 
(a) partial contribution from each band and (b) total field 
distribution. The k = \/8'Ke^Vp^NFcr /c^ /£,cr that comes into 

14 for MgB2). 



the Maxwell equation is taken to be 10 (k : 



D. Current density 

In figure [T^ the current densities j(r) in the two bands 
as a function of distance r from the vortex centre are 
shown for ^Tr/^a = 3, po = 0.5, = ria, A ~ —0.1, 
and Tcji^a — 10, for various temperatures. The current 
density contribution from the tt band can be substantial, 
and even dominating as demonstrated in this figure, for 
larger po and ^tt/Co-- As temperature is lowered, the a- 
band current density exhibits the KP effect and is more 
confined around the vortex centre. For stronger coupling 
of the two bands, the KP effect becomes manifest also in 
the TT band. 

Deviation of the magnetic field distribution from the 
uniform field, SH{r) ~ H{r) — Hq, can be obtained from 
j(r) by integrating the Maxwell equation, where a pa- 
rameter K = y^Sne^Vp^Npa / / comes in. For MgB2, 
K « 14. The SH{r) as a function of r is illustrated in 
figuremfor = 3, = n„ A = -0.1, T/T, = 0.1, 

Tc/S.a — 10, and K = 10, for po — 0.3 and 0.5. In fig- 
ure I13r a) partial contribution from each band is shown, 
while [T^b) presents the total field distribution. It can 
be seen that when coupling of the two bands is increased, 
field fluctuations in the n band (hence those in the total 
field) are enhanced significantly, while the a band con- 
tribution hardly changes. 




FIG. 14: Upper critical field as a function of temperature T 
for Po = 0.3, n-n = n<T, A = ±0.1, for (a) ^n/^rr ~ 3 and 
(b) 5. The result for a single clean band is also shown. The 
hc/\e\(l 13 T for MgBa. 



E. Phase diagram 

We show in figure [H] the upper critical field Hc2 (T) 
for T near T, for po = 0.3, n^r — ria-, and A — ±0.1, 
for (a) Ctt/Co- = 3 and (b) 5, along with the result for 
a single clean band. Although Hc2 is sensitive to the 
amount of impurities in the sample, in MgB2 single crys- 
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tals, Hc2 for field parallel to the c axis has been mea- 
sured to be HrojT = TJ2\k, 2 Tesla in several exper- 
iments [13, [U, M [13, 0, [si, [EE m [si, [si, [el . This 

value roughly corresponds to A somewhere between -0.1 
and 0.1 for ^tt/'^o- — Pa — 0.3, and — n„ in our 
model. When the a band is coupled with the tt band, 
Hc2 is substantially reduced from the single-band value, 
and decreases further when the Coulomb repulsion dom- 
inates or the diffusivity increases in the tt band. This 
is consistent with stronger suppression of the order pa- 
rameter and enlargement of the vortex core area in the a 
band. Furthermore, for relatively large (,Tr/(,a, Hc2{T) de- 
velops an upward curvature near Tc (figure [TWb)). This 
is interesting in light of the theoretical work based on the 
Eliashberg theory by Mansor and Carbotte, who studied 
the effects of Fermi velocity anisotropy and impurities 
on Hc2 Their prediction for MgB2 is that, for field 
in the c direction (where the Fermi velocities in the two 
bands are assumed to be the same), Hc2(T) exhibits an 
upward curvature near when the tt band is clean and 
the a band is dirty, while it has a quasilinear T depen- 
dence for the dirty tt and the clean u band (as observed 
in several experiments; see references in ref. (6]|). Such 
an upward curvature in Hc2 (|| c) has been observed in 
some experiments on MgB2 single crystals [H, [gl, [g^l • 

IV. CONCLUSION 



critical temperature and field (where the order parameter 
at the vortex unit cell boundary vanishes) are reduced, 
and increased diffusivity in the dirty band can result in 
an upward curvature of the upper critical field near the 
transition temperature. The Kramer-Pesch effect arising 
from thermal depopulation of higher-energy core bound 
states tends to disappear as field becomes stronger and 
overlap of quasiparticle states of neighbouring vortices in- 
creases. The zero-bias LDOS in the diffusive band has a 
decay length much larger than that in the ballistic band, 
indicating substantial intervortex transfer of quasiparti- 
cles, and the half-decay length increases significantly as 
the diffusivity increases. Coupling with the "weak" diffu- 
sive band leads to smearing of energy bands of the vortex 
lattice and van Hove singularities at the band edges in 
the vortex-core LDOS in the "strong" ballistic band. Fur- 
thermore, bound states tend to appear at the gap edge in 
the ballistic band, in addition to the well-known Caroli- 
de Gennes-Matricon bound states. These effects are en- 
hanced for increased coupling, diffusivity, and Coulomb 
repulsion. Finally, the current density contribution in 
the vortex core and resulting field fluctuations in the dif- 
fusive band can be substantial, or even dominating, and 
exhibit the Kramer-Pesch effect manifest in the current 
density when coupling with the ballistic band is strong. 
We find the above intriguing features in the quasiparticle 
spectra for parameter values appropriate for MgB2 . 



In conclusion, we have studied the effects of induced 
superconductivity and impurities on the electronic struc- 
ture in the vortex lattice of a two-band superconductor, 
in which a "weak" 3D diffusive band and a "strong" 2D 
ballistic band are hybridized through the pairing interac- 
tion. We have found that the Coulomb repulsion and the 
diffusivity in the "weak" dirty band enhance suppression 
of the order parameter and enlargement of the vortex core 
by magnetic field in the "strong" clean band. As a result, 
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